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PROVISIONAL APPLICATION FOR PATENT 



Solution Additives for the Attenuation of 
Protein Aggregation 



Background of the Invention 

The process of protein folding is complex, and a complete understanding of it is one of 
the grand challenges facing contemporary biochemists. This complexity arises in part from the 
fact that a nascent protein may not simply fold into its native state under the influence of only the 
primary solvent (water), but may also interact with other molecules in solution. The effects of 
other molecules may be favorable for folding, as is the case for molecules like folding 
chaperones, or unfavorable, as is the case for other partially-unfolded protein molecules. 

One of the primary driving forces in protein folding is the burial of exposed hydrophobic 
residues. Dill, K. A. Biochemistry 1990, 29, 7133-7155. Aggregation results if the hydrophobic 
collapse should occur in an intermolecular instead of an intramolecular fashion. Because 
aggregation occurs as a parallel reaction to proper folding, there is kinetic competition between 
the two pathways. Orsini, G.; Goldberg, VI. E. J. Biol. Chem. 1978, 253, 3453-3458; 
Zettlmeissl, G.; Rudolph; R.; Jaenicke. R. Biochemistry 1979, 18, 5567-5571; Kiefllaber, T.; 
Rudolph; R.; Kohler, H.-H.; Buchner, J. Bio/Technology 1991, 9, 825-829; Hevehan, D. L.; 
Clark, E. D. B. Biotechnol. Bioeng. 1997, 54, 221-230. 

Aggregation of misfolded proteins is a s ignificant p roblem b oth / « vivo and in vitro. 
Aggregation has been implicated in human diseases, such as Huntington's, Alzheimer's, and 
Parkinson's Diseases. Taylor, J. P.; Hardy, J.; Fischbeck; K. H. Science 2002, 296, 1991-1995. In 
applied biotechnology, aggregation is a significant side reaction of protein refolding, which is an 
important step in the production of many recombinant proteins. De Bemandez Clark, E.; 
Schwarz, E.; Rudolph, R. Methods Enzymol. 1999, 309, 217-236. 

Both nature and man have developed strategies to combat aggregation. Chaperonins, 
such as the GroEL/GroES system, surround and isolate partially-folded proteins in the bulk 
cytosol so they can continue to fold without aggregating. Hartl, F. U.; Hayer-Hartl, M. Science 
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2003, 295, 1852-1858. Similarly, additives to deter aggregation are often included in protein 
refolding buffers and other in vitro applications, such as pharmaceutical formulations. Wang, 
W.Int. J.Pharm. 1999, 185, 129-188. 



Presently disclosed are classes of additives that, when added to protein solutions, 
attenuate the rate of aggregation. The members of the classes have two key, well-defined 
properties that result in their ability to slow aggregation. The present invention also 
recognizes that there are many other such molecules that have never been synthesized and 
have never been used to stabilize proteins. Some of them exemplify the two properties above 
to a greater extent and are, therefore, superior solution additives for the attenuation of protein 
aggregation. 

In one embodiment the present invention relates to a compound comprising a non- 
protein-binding moiety (NPBM) and at least one protein binding group (PBG). In a further 
embodiment, the NPBM is a polyol, sugar, amino acid, or dendrimer moiety. In a further 
embodiment, the polyol moiety is a sorbitol or mannitol moiety. In a further embodiment, the 
sugar moiety is a glucose, sucrose, or trehalose moiety. In a fiirther embodiment, the amino acid 
moiety is an arginine betaine, proline, or ectoine moiety. In a further embodiment, the dendrimer 
moiety is based on benzene, pentaerythritol, P(CH20H)3, or TRIS. 

In a further embodiment, the PBG is a urea, guanidinium ion, detergent, amino acid, 
denaturant, surfactant, polysorbate, polaxamer, citrate, chaotrope, or acetate group. In a further 
embodiment, the PBG is a guanidinium ion. In a further embodiment, the PBG is sodium 
dodecyl sulfate. 

In another embodiment, the present invention relates to a compound of formula I: 



Summary of the Invention 




I 



wherein: 
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R is an electron pair, H, alkyl, aryl, heteroaryl, aralkyl, heteroaralkyl, or an alkali metal; 

R' is H, alkyl, aryl, heteroaryl, aralkyl, heteroaralkyl, or R"H2N; 

R" is an electron pair, H, alkyl, aryl, heteroaryl, aralkyl, or heteroaralkyl; 

W is O, NH2*(halogen)", or S; and 

nis 1,2, or 4-100. 

In a further embodiment, the present invention relates to a compound of formula I and the 
attendant d efmitions, w herein R i s an e lectron p air. hi a f urther e mbodiment, R ' i s H . I n a 
further embodiment, R' is R"H2N. In a further embodiment, R' is H3N*. In a further 
embodiment, W is NH2*Cr. In a further embodiment, nisi. In a fiirther embodiment, n is 2. In 
a further embodiment, n is 4. In a further embodiment, n is 5. In a further embodiment, n is 6. 
In a further embodiment, R is an electron pair, R' is Hibt, W is NH2^Cr, and nisi. In a further 
embodiment, R is an electron pair, R' is HsIsT, W is NH2^Cr, and n is 2. In a further 
embodiment, R is an electron pair, R' is HsN"^, W is NH2*Cr, and n is 4. In a further 
embodiment, R is an electron pair, R' is H3N*, W is NH2*Cr, and n is 5. In a further 
embodiment, R is an electron pair, R' is HaN^, W is NH2'^Cr, and n is 6. In a further 
embodiment, R is an electron pair, R' is HsN^, W is O, and n is 1. In a fiirther embodiment, R is 
an electron pair, R' is HaN^, W is O, and n is 2. In a further embodiment, R' is Kibt, W is O, 
and n is 4. In a further embodiment, R is an electron pair, R' is HaN^, W is O, and n is 5. In a 
further embodiment, R is an electron pair, R' is HjN^, W is O, and n is 6. In a further 
embodiment, R is an electron pair, R' is H, W is NH2*Cr, and n is 1 . In a further embodiment, R 
is an electron pair, R' is H, W is NH2^Cr, and n is 2. In a further embodiment, R is an electron 
pair, R' is H*, W is NH2*Cr, and n is 4. In a further embodiment, R is an electron pair, R' is H, 
W is NH2^Cr, and n is 5. In a further embodiment, R is an electron pair, R' is H, W is NH2^Cr, 
and n is 6. In a further embodiment, R is an electron pair, R' is H, W is O, and n is 1. In a 
further embodiment, R is an electron pair, R' is H, W is O, and n is 2. In a further embodiment, 
R is an electron pair, R' is H, W is O, and n is 4. In a further embodiment, R is an electron pair, 
R' is H, W is O, and n is 5. In a further embodiment, R is an electron pair, R' is H, W is O, and 
n is 6. 

In another embodiment, the present invention relates to one of the following compounds: 
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wherein, independently for each occurrence, 
RisHorCH2Y; 

R' is H, a sugar radical, or CH2 Y; 

n is an integer from 1 to 100, inclusive; 

a is 1, 2, or 3; 

X is C(CH2Y)3; and 

Y is a protein binding group, 

wherein at least one Y is present in all compounds. 

In a further embodiment, Y is a guanidinium ion. 



-4- 



In another embodiment, the present invention relates to a method of screening 
compounds for the property of inhibiting protein aggregation in solution, comprising: 

a) computing a set of parameters utilizing molecular modeling based on compounds 
known to have the property of inhibiting protein aggregation; 

b) applying those parameters to other compounds; and 

c) choosing the compounds that meet the criteria of those parameters. 

In another embodiment, the present invention relates to a method of preparing new 
compounds having the property of protein aggregation inhibition in solution, comprising: 

a) computing a set of parameters utilizing molecular modeling based on compounds 
known to have the property of inhibiting protein aggregation; 

b) designing compounds based on those parameters; and 

c) synthesizing the compounds. 

In another embodiment, the present invention relates to a method of classifying 
compounds a s e ither i nhibitory o f p rotein a ggregation i n s olution o r n ot i nhibitory o f p rotein 
aggregation in solution, comprising: 

a) determining the phase space trajectories of the protein, solvent, and additive using 
molecular dynamics; 

b) calculating the distance, r, between the center of mass for both the solvent molecule 
and additive molecule to the protein's van der Waals surface; 

c) determining the minimum distance, r*, at which no significant differences between the 
local (r = r*) and bulk density are observed; 

d) determining which molecules lie within the distance, r*, from the protein surface and 
classifying these molecules as the local domain; 

e) determining which molecules lie outside the distance, r*, from the protein surface and 
classifying these molecules as the bulk domain; 

f) determining the instantaneous preferential binding coefficient, rxp(t), using the 
following formula: 
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rxp(t) = n"x - n'x (n"w / n'w) 



wherein: 

n"x = the number of additive molecules in the bulk domain; 
n'x = the number of additive molecules in the local domain; 
n"w = the number of solvent molecules in the bulk domain; and 
n'w = the number of solvent molecules in the local domain; and 

g) calculating the preferential binding coefficient, Fxp, as the time average of each of the 
values in step f) using the following formula: 



These embodiments of the present invention, other embodiments, and their features and 
characteristics, will be apparent from the description, drawings and claims that follow. 



Figure 1 depicts a simplified dimerization reaction-coordinate diagram for the reaction 
U + U — * A2 (equation 2). The dotted line is the reaction coordinate in water and the solid 
line is the reaction coordinate in the presence of an additive having the two anti-aggregation 
properties d iscussed. Protein molecules are represented by black coils and the additive by 
dark grey circles. The energy difference between the reactants (U + U) and the t ransition 
state determines the rate of the reaction. In the A2 state, the region between the protein 
molecules (light grey oval) is preferentially hydrated because water can enter this region but 
the additive c annot. This preferential h ydration i ncreases the free e nergy o f t he t ransition 
state, increases the energy barrier for the reaction, and slows the reaction rate. 

Figure 2 depicts arginine derivatives with shorter (left) and longer (right) methylene 
linkers between their amino acid backbone and guanidino functional groups. 

Figure 3 depicts molecules that will be preferentially-oriented at the protein-solvent 
interface. Molecule ( a) i s a d erivative o f g lucose ( stabilizer) 1 inked to a d imethyl-guanidino 




0 



Brief Description of the Figures 
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(destabilizer) moiety. Molecule (b) is a polyol (stabilizer) with a guanidino group 
(destabilizer) attached to one end. 

Figure 4 depicts the physical interpretation of the preferential binding coefficient. 
Interactions of solvent molecules with the protein at the protein-solvent interface generally 
induce solvent concentration differences in the local (II) and bulk (I) domains, fxp is the 
thermodynamic measure of the number of additive molecules boimd to the protein, or in other 
words, the excess number of additive molecules in the vicinity of the protein versus the number 
of additive molecules in an equivalent volume of bulk solution. 

Figure 5 depicts a simulation cell containing RNase Tl (center spheres) solvated by 
water (thin lines) and urea (spheres). 

Figure 6 depicts radial distribution functions of water, urea, and glycerol shown for 
simulations of RNase Tl in glycerol and urea solutions (left) and RNase A in a glycerol solution 
(right). In the left-hand figure, the difference between the two gw(r) functions is not visible at 
this scale. 

Figure 7 depicts apparent preferential binding coefficient as a function of the cutoff 
distance between the local and bulk domains for simulations of RNase Tl in glycerol and urea 
solution. 

Figure 8 depicts rxp(/) probability density fimction. A wide range of values of rxp(0 are 
sampled as water and cosolvent molecules diffuse between the local and bulk domains. 

Figure 9 depicts the correlation of solvent-accessible area and the number of water 
molecules in the local domain of constituent groups. Each point represents a constituent group 
of either a type of amino acid side chain or the protein backbone in one of the three simulations 
shown in Table 2. The solvent accessible area of a constituent group and the number of water 
molecules in the local domain of the solvent near the group (riwi) are correlated. 

Figure 10 depicts the binding behavior of glycerol and water with the 15 serine residues 
in RNase Tl as shown in a plot of the number of glycerol molecules in the local domain of each 
serine residue versus the number of water molecules in the same volume. The labels are the one- 
letter codes for each amino acid side chain, and "B" is the protein backbone. The line represents 
the bulk glycerol composition. Ser 17, 35, and 72 have positive preferential binding coefficients. 
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Ser 63 has a negative preferential binding coefficient, and the remaining 1 1 serine residues have 
essentially zero values for their preferential binding coefficients. 

Figure 11 depicts the local binding behavior of urea and water with the amino acid 
backbone and side chains in RNase Tl. The labels are the one-letter codes for the amino acid 
side chains, and "B" is the protein backbone. The line denotes the bulk urea concentration. In 
addition to the protein backbone and Ser, the hydrophobic amino acids Cys, Gly, Leu, Phe, Pro, 
Tyr, and Val all preferentially bind urea, while the hydrophilic Asp preferentially binds water. 

Figure 12 depicts the group preferential binding coefficients for glycerol with the amino 
acid backbone and side chains in RNase Tl. The labels are the one-letter codes for the amino 
acid side chains, and "B" is the protein backbone. The line denotes the bulk glycerol 
concentration. Tyr and Gly preferentially bind glycerol; Asp and Glu preferentially bind water; 
and the binding coefficients of the other groups are not statistically different fi-om zero. 

Figure 13 depicts the local binding behavior of glycerol with the amino acid backbone 
and side chains in RNase A. The labels are the one-letter codes for the amino acid side chains, 
and "B" is the protein backbone. The line denotes the bulk glycerol concentration. All of the 
constituent groups in RNase A either preferentially bind water or are neutral. 

Detailed Description of the Invention 

Definitions 

For convenience, before fiirther description of the present invention, certain terms 
employed in the specification, examples and appended claims are collected here. These 
definitions s hould b e r ead i n 1 ight o f t he r emainder o f t he d isclosure a nd u nderstood a s b y a 
person of skill in the art. Unless defined otherwise, all technical and scientific terms used herein 
have the same meaning as commonly understood by a person of ordinary skill in the art. 

The articles "a" and "an" are used herein to refer to one or to more than one (i.e., to at 
least one) of the grammatical object of the article. By way of example, "an element" means one 
element or more than one element. 

The terms "comprise" and "comprising" are used in the inclusive, open sense, meaning 
that additional elements may be included. 
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The term "including" is used to mean "including but not limited to". "Including" and 
"including but not limited to" are used interchangeably. 

The term "additive" as used herein refers to any component other than the subject protein 
and the main solvent. Non-limiting examples of additives include small molecules, cosolvents, 
buffer salts, and stabilizers. 

The term "dendrimer" is used to mean a broad class of polymers constructed via stepwise 
polymerization from a central "core unit," one or more "branching units," and several "surface 
units." The review of Matthews (1998) provides an overview of dendrimers including 
compositions and synthetic routes. Core units may include (but are not limited to) carbon, 
nitrogen, phosphorous, benzene, and porphyrins. A non-extensive collection of 17 specific 
chemistries that are used to create branching units are sunmiarized in Table 2 of Matthews 
(1998). 

The term "TRIS" is art-recognized and refers to tris(hydroxymethyl)aminomethane. 

The term "aliphatic" is an art-recognized term and includes linear, branched, and cyclic 
alkanes, alkenes, or alkynes. In certain embodiments, aliphatic groups in the present invention 
are linear or branched and have from 1 to about 20 carbon atoms. 

The term "alkyl" is art-recognized, and includes saturated aliphatic groups, including 
straight-chain alkyl groups, branched-chain alkyl groups, cycloalkyl (alicyclic) groups, alkyl 
substituted cycloalkyl groups, and cycloalkyl substituted alkyl groups. In certain embodiments, a 
straight chain or branched chain alkyl has about 30 or fewer carbon atoms in its backbone (e.g., 
C1-C30 for straight chain, C3-C30 for branched chain), and alternatively, about 20 or fewer. 
Likewise, cycloalkyls have from about 3 to about 10 carbon atoms in their ring structure, and 
alternatively about 5, 6 or 7 carbons in the ring structure. 

Unless the number of carbons is otherwise specified, "lower alkyl" refers to an alkyl 
group, as defined above, but having from one to ten carbons, alternatively from one to about six 
carbon atoms in its backbone structure. Likewise, "lower alkenyl" and "lower alkynyl" have 
similar chain lengths. 

The term "aralkyl" is art-recognized, and includes alkyl groups substituted with an aryl 
group (e.g., an aromatic or heteroaromatic group). 
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The terms "alkenyl" and "alkynyl" are art-recognized, and include unsaturated aliphatic 
groups analogous in length and possible substitution to the alkyls described above, but that 
contain at least one double or triple bond respectively. 

The term "heteroatom" is art-recognized, and includes an atom of any element other than 
carbon or hydrogen. Illustrative heteroatoms include boron, nitrogen, oxygen, phosphorus, sulfur 
and selenium, and alternatively oxygen, nitrogen or sulfur. 

The term "aryl" is art-recognized, and includes 5-, 6- and 7-membered single-ring 
aromatic groups that may include from zero to four heteroatoms, for example, benzene, 
naphthalene, anthracene, pyrene, pyrrole, furan, thiophene, imidazole, oxazole, thiazole, triazole, 
pyrazole, pyridine, pyrazine, pyridazine and pyrimidine, and the like. Those aryl groups having 
heteroatoms in the ring structure may also be referred to as "heteroaryl" or "heteroaromatics." 
The aromatic ring may be substituted at one or more ring positions with such substituents as 
described above, for example, halogen, azide, alkyl, aralkyl, alkenyl, alkynyl, cycloalkyl, 
hydroxyl, alkoxyl, amino, nitro, sulfhydryl, imino, amido, phosphonate, phosphinate, carbonyl, 
carboxyl, silyl, ether, alkylthio, sulfonyl, sulfonamido, ketone, aldehyde, ester, heterocyclyl, 
aromatic or heteroaromatic moieties, -CF3, -CN, or the like. The term "aryl" also includes 
polycyclic ring systems having two or more cyclic rings in which two or more carbons are 
common to two adjoining rings (the rings are "fused rings") wherein at least one of the rings is 
aromatic, e.g., the other cyclic rings may be cycloalkyls, cycloalkenyls, cycloalkynyls, aryls 
and/or heterocyclyls. 

The terms ortho. meta and para are art-recognized and apply to 1,2-, 1,3- and 1,4- 
disubstituted benzenes, respectively. For example, the names 1,2-dimethylbenzene and ortho- 
dimethylbenzene are synonymous. 

The terms "heterocyclyl" and "heterocyclic group" are art-recognized, and include 3- to 
about 1 0-membered ring structures, such as 3- to about 7-membered rings, whose ring structures 
include one to four heteroatoms. Heterocycles may also be polycycles. Heterocyclyl groups 
include, for example, thiophene, thianthrene, furan, pyran, isobenzofuran, chromene, xanthene, 
phenoxathiin, pyrrole, imidazole, pyrazole, isothiazole, isoxazole, pyridine, pyrazine, 
pyrimidine, pyridazine, indolizine, isoindole, indole, indazole, purine, quinolizine, isoquinoline, 
quinoline, phthalazine, naphthyridine, quinoxahne, quinazoline, cinnoline, pteridine, carbazole, 
carboline, phenanthridine, acridine, pyrimidine, phenanthroline, phenazine, phenarsazine. 
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phenothiazine, furazan, phenoxazine, pyrrolidine, oxolane, thiolane, oxazole, piperidine, 
piperazine, morpholine, lactones, lactams such as azetidinones and pyrrolidinones, sultams, 
sultones, and the like. The heterocyclic ring may be substituted at one or more positions with 
such substituents as described above, as for example, halogen, alkyl, aralkyl, alkenyl, alkynyl, 
cycloalkyl, hydroxyl, amino, nitro, sulfhydryl, imino, amido, phosphonate, phosphinate, 
carbonyl, carboxyl, silyl, ether, alkylthio, sulfonyl, ketone, aldehyde, ester, a heterocyclyl, an 
aromatic or heteroaromatic moiety, -CF3, -CN, or the like. 

The terms "polycyclyl" and "polycyclic group" are art-recognized, and include structures 
with two or more rings (e.g., cycloalkyls, cycloalkenyls, cycloalkynyls, aryls and/or 
heterocyclyls) in which two or more carbons are common to two adjoining rings, e.g., the rings 
are "fused rings". Rings that are joined through non-adjacent atoms, e.g., three or more atoms are 
common to both rings, are termed "bridged" rings. Each of the rings of the polycycle may be 
substituted with such substituents as described above, as for example, halogen, alkyl, aralkyl, 
alkenyl, alkynyl, cycloalkyl, hydroxyl, amino, nitro, sulfhydryl, imino, amido, phosphonate, 
phosphinate, carbonyl, carboxyl, silyl, ether, alkylthio, sulfonyl, ketone, aldehyde, ester, a 
heterocyclyl, an aromatic or heteroaromatic moiety, -CF3, -CN, or the like. 

The term "carbocycle" is art-recognized and includes an aromatic or non-aromatic ring in 
which each atom of the ring is carbon. The flowing art-recognized terms have the following 
meanings: "nitro" means -NO2; the term "halogen" designates -F, -CI, -Br or -I; the term 
"sulfhydryl" means -SH; the term "hydroxyl" means -OH; and the term "sulfonyl" means -SO2'. 

The terms "amine" and "amino" are art-recognized and include both unsubstituted and 
substituted amines, e.g., a moiety that may be represented by the general formulas: 

R50 

R53 

\ I 

R52 

wherein R50, R51 and R52 each independently represent a hydrogen, an alkyl, an alkenyl, - 
(CH2)m-R61, or R50 and R51, taken together with the N atom to which they are attached 
complete a heterocycle having from 4 to 8 atoms in the ring structure; R61 represents an aryl, a 
cycloalkyl, a cycloalkenyl, a heterocycle or a polycycle; and m is zero or an integer in the range 
of 1 to 8. In certain embodiments, only one of R50 or R51 may be a carbonyl, e.g., R50, R51 and 
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the nitrogen together do not form an imide. In other embodiments, R50 and R5 1 (and optionally 
R52) each independently represent a hydrogen, an alkyl, an alkenyl, or -(CH2)ni-R61. Thus, the 
term "alkylamine" includes an amine group, as defined above, having a substituted or 

unsubstituted alkyl attached thereto, i.e., at least one of R50 and R51 is an alkyl group. 

The term "acylamino" is art-recognized and includes a moiety that may be represented by 
the general formula: 

0 



-N- 



R50 



-R54 



wherein R50 is as defined above, and R54 represents a hydrogen, an alkyl, an alkenyl or - 
(CH2)m-R61, where m and R61 are as defined above. 

The term "amido" is art-recognized as an amino-substituted carbonyl and includes a 
moiety that may be represented by the general formula: 

O 

R51 



R50 

wherein R50 and R51 are as defined above. Certain embodiments of the amide in the present 
invention will not include imides which may be unstable. 

The term "alkylthio" is art-recognized and includes an alkyl group, as defined above, 
having a sulfur radical attached thereto. In certain embodiments, the "alkylthio" moiety is 
represented by one of -S-alkyl, -S-alkenyl, -S-alkynyl, and -S-(CH2)m-R61, wherein m and R61 
are defined above. Representative alkylthio groups include methylthio, ethyl thio, and the like. 

The term "carbonyl" is art-recognized and includes such moieties as may be represented 
by the general formulas: 



O 



O 



R55 



X50 X50 R56 

wherein X50 is a bond or represents an oxygen or a sulfiu", and R55 represents a hydrogen, an 
alkyl, an alkenyl, -(CH2)m-R61or a pharmaceutically acceptable salt, R56 represents a hydrogen, 
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an a Ikyl, an a Ikenyl o r -(CH2)m-R6 1 , w here m and R 6 1 a re d efined above. W here X 50 i s a n 
oxygen and R55 or R56 is not hydrogen, the formula represents an "ester". Where X50 is an 
oxygen, and R55 is as defined above, the moiety is referred to herein as a carboxyl group, and 
particularly when R55 is a hydrogen, the formula represents a "carboxylic acid". Where X50 is 
an oxygen, and R56 is hydrogen, the formula represents a "formate". In general, where the 
oxygen atom of the above formula is replaced by sulfur, the formula represents a "thiocarbonyl" 
group. Where X50 is a sulfur and R55 or R56 is not hydrogen, the formula represents a 
"thioester." Where X50 is a sulfur and R55 is hydrogen, the formula represents a "thiocarboxylic 
acid." Where X50 is a sulfur and R56 is hydrogen, the formula represents a "thioformate." On 
the other hand, where X50 is a bond, and R55 is not hydrogen, the above formula represents a 
"ketone" group. Where X50 is a bond, and R55 is hydrogen, the above formula represents an 
"aldehyde" group. 

The terms "alkoxyl" or "alkoxy" are art-recognized and include an alkyl group, as 
defined above, having an oxygen radical attached thereto. Representative alkoxyl groups include 
methoxy, ethoxy, propyloxy, tert-butoxy and the like. An "ether" is two hydrocarbons covalently 
linked by an oxygen. Accordingly, the substituent of an alkyl that renders that alkyl an ether is or 

resembles an alkoxyl, such as may be represented by one of -0-alkyl, -0-alkenyl, -0-alkynyl, - 
0-(CH2)m-R61, where m and R61 are described above. 

The term "sulfonate" is art-recognized and includes a moiety that may be represented by 
the general formula: 

O 

S OR57 

O 

in which R57 is an electron pair, hydrogen, alkyl, cycloalkyl, or aryl. 

The term "sulfate" is art-recognized and includes a moiety that may be represented by the 
general formula: 

O 

I 

O S OR57 

O 
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in which R57 is as defined above. 

The term "sulfonamido" is art-recognized and includes a moiety that may be represented 
by the general formula: 

O 



-N- 



-OR56 



R50 O 

in which R50 and R56 are as defined above. 

The term "sulfamoyl" is art-recognized and includes a moiety that may be represented by 
the general formula: 

O 



-N, 



/ 
\ 



R50 



O 



R51 



in which R50 and R5 1 are as defined above. 

The term "sulfonyl" is art-recognized and includes a moiety that may be represented by 
the general formula: 



1 


1 


1 


5 ] 

1 


1 


II 

0 



R58 



in which R58 is one of the following: hydrogen, alkyl, alkenyl, alkynyl, cycloalkyl, heterocyclyl, 
aryl or heteroaryl. 

The term "sulfoxido" is art-recognized and includes a moiety that may be represented by 
the general formula: 



\ 



P 



R58 

in which R58 is defined above. 

The term "phosphoramidite" is art-recognized and includes moieties represented by the 
general formulas: 
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o 



Q51 — p 



-0- 



N 
/ \ 
R50 R51 



0 



Q51 — p— OR59 



R50 R51 



wherein Q51, R50, R51 and R59 are as defined above. 

The term "phosphonamidite" is art-recognized and includes moieties represented by the 
general formulas: 

R60 R60 



Q51_p 



-0- 



R50 R51 



Q51 — p— OR59 



R50 R51 



wherein Q51, R50, R51 and R59 are as defined above, and R60 represents a lower alkyl or an 
aryl. 

Analogous substitutions may be made to alkenyl and alkynyl groups to produce, for 
example, aminoalkenyls, aminoalkynyls, amidoalkenyls, amidoalkynyls, iminoalkenyls, 
iminoalkynyls, thioalkenyls, thioalkynyls, carbonyl-substituted alkenyls or alkynyls. 

The definition of each expression, e.g. alkyl, m, n, etc., when it occurs more than once in 
any structure, is intended to be independent of its definition elsewhere in the same structure 
unless otherwise indicated expressly or by the context. 

For purposes of this invention, the chemical elements are identified in accordance with 
the Periodic Table of the Elements, CAS version. Handbook of Chemistry and Physics, 67th Ed., 
1986-87, inside cover. 

Overview 

Proteins are widely used in medical and industrial applications. One of the major 
difficulties encountered in these applications is that proteins are prone to degradation by a 
variety of routes, the most common of which is aggregation. Aggregated protein generally 
does not have the same functionality as normal, native protein. For example, in 
pharmaceutical applications, the consequences of administering aggregated drug to a patient 
can be severe because aggregates can be cytotoxic; and they generally induce an immune 
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response. Bucciatini, M.; Giannoni, E.; Chiti, F.; Baroni, F.; Formigh, L.; Zurdo, J.; Taddei, 
N.; Ramponi, G.; Dobson, C. M.; Stefani, M. Nature 2002, 416, 507-51 1; Braun, A.; Kwee, L.; 
Labow, M. A.; Alsenz, J. Pharm. Res. 1997, 14, 1472-1478. Due to these and other negative 
effects, protein solutions often contain one or more additives designed to deter aggregation. 
Wang, W. Int. J. Pharm. 1999, 185,129-188. 

As a protein folds, there is kinetic competition between the proper folding pathway 
and side reactions, such as aggregation. This is depicted schematically in equations 1 and 2: 

U-^N (1) 

U + U A2 (2) 

where U represents an unfolded protein; N represents a folded, native protein; and A2 
represents a small aggregate species. 

Alternatively, if the protein is initially in its native state, such as in a pharmaceutical 
formulation, aggregation proceeds through formation of a partially-unfolded intermediate, I, 
which can aggregate in a sense analogous to an unfolded protein: 

N ^ I (3) 

I + 1 ^ A2 (4) 

For industrial and medical applications, it is desirable to eliminate or minimize the 
formation of protein aggregates. In protein folding or refolding processes, decreasing the rate 
of aggregation results in a higher yield of active, properly-folded protein. In pharmaceutical 
formulations, decreasing the rate of aggregation causes more drug to remain in its active form 
and eliminates the possibly dangerous side effects of administering aggregated protein to the 
patient. To minimize aggregation, various conditions, such as temperature, pH, and the type 
and amount of buffer additives, are screened experimentally to identify an optimum set of 
conditions. 

Additives that attenuate protein aggregation have historically been identified by 
heuristic, experimental screening methods. Identification of specific additive-protein 
combinations that result in decreased aggregation via these techniques has led to little 
fundamental understanding about what properties of a molecule lead to the ability to 
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suppress protein aggregation. This gap in understanding has prevented development of 
rational strategies to prevent protein aggregation. 

Through the mechanistic understanding summarized presently, two fundamental 
properties of a good anti -aggregation additive have been identified. This discovery allows 
additives to be selected based on their relative ranking in terms of these two properties, thus 
narrowing experimental testing to molecules likely to have optimal performance. It also 
enables molecules to be classified based on whether they may have the ability to attenuate 
aggregation. The rational, mechanistic classification schemes of the present invention will 
allow entire classes of protein-aggregation-attenuating additives and formulations to be 
identified. 

Additionally, a quantitative method based on molecular dynamics simulations using all 
atom potential models has been developed and validated for calculating preferential binding 
coefficients. The present invention is not a derivative of thermodynamic integration or 
thermodynamic perturbation methods and requires only a single trajectory to compute the 
transfer free energy of a protein into a weak-binding additive system. The results match 
experimental data well for glycerol and urea solutions, covering a range of positive and negative 
binding behavior. The present invention also augments experimentally-observable, macroscopic 
thermodynamics with the mechanistic insight provided by a molecular-level, statistical 
mechanical model. 

Variations in the radial distribution functions with distance for each additive are evident 
up to about 6A, i.e., roughly two solvation shells of water, away from the protein. Glycerol is 
not totally excluded from close contact with the protein, but glycerol is less likely than urea to be 
found in such a position. The radial distribution functions of water and additives are sufficient to 
calculate preferential binding coefficients by integrating over a suitable solvent volume. 

The binding behavior of the amino acid side chains in RNase Tl qualitatively follow a 
hydrophilic series, with more hydrophilic amino acids in the protein tending to have a higher 
concentration of water in their vicinity. The constituent group binding behavior differs between 

the groups in RNase A to those in RNase Tl. Development of a group contribution method at the 
amino acid level for estimating binding coefficients or transfer fr-ee energies of whole proteins is 
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complicated by the wide range of coordination behaviors observed for single types of amino 
acids in different environments on the protein surface. 

In the pharmaceutical industry, many protein drugs are synthesized in bacterial hosts, 
such as E. coli, in the form of solid, partially-aggregated precipitates called inclusion bodies. 
These inclusion bodies must be unfolded and solubilized, and then refolded to form active 
protein. During refolding, proteins are especially susceptible to aggregation, and additives 
must be used to minimize aggregation and increase the yield of biologically-active protein. 
The compounds 0 f the present invention are ideal for use in these circumstances because 
they will slow the rate of aggregation and therefore increase the yield of active protein. 
Likewise, when pharmaceutically-active proteins are formulated in aqueous solution, 
additives are used to prevent aggregation during storage, thereby increasing its shelf-life. 
The compounds of the present invention are also useful in preventing aggregation in these 
circumstances. Additional applications can be envisioned by those of ordinary skill in the art of 
protein stabilization. The above applications are meant to be only exemplary and not limiting in 
any way. 

Select Preferred Embodiments 

In a preferred embodiment, the present invention relates to a method of suppressing or 
preventing aggregation of a protein in solution, comprising the step of combining in a solution 
a compound of the present invention and a protein. In certain embodiments, the protein is a 
recombinant protein. In certain embodiments, the protein is a recombinant antibody. In 
certain embodiments, the protein is a recombinant human antibody. In certain embodiments, 
the protein is a recombinant mammalian protein. In certain embodiments, the protein is a 
recombinant human protein. In certain embodiments, the protein is recombinant human 
insulin, recombinant human erythropoietin or a recombinant human interferon. In certain 
embodiments, the so lution i s a n a queous so lution. I n c ertain e mbodiments, the p rotein is a 
recombinant protein; and the solution is an aqueous solution. In certain embodiments, the 
protein is a recombinant human antibody; and the solution is an aqueous solution. In certain 
embodiments, the protein is a recombinant human protein; and the solution is an aqueous 
solution. 
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In a preferred embodiment, the present invention relates to a method of suppressing or 
preventing aggregation of a protein in solution, comprising the step of combining in a solution 
a compound of the present invention and a protein, hi certain embodiments, the protein is a 
recombinant protein. In certain embodiments, the protein is a recombinant antibody. In 
certain embodiments, the protein is a recombinant human antibody. In certain embodiments, 
the protein is a recombinant mammalian protein. In certain embodiments, the protein is a 
recombinant human protein. In certain embodiments, the protein is recombinant human 
insulin, recombinant human erythropoietin or a recombinant human interferon. In certain 
embodiments, the so lution i s a n a queous so lution. I n c ertain e mbodiments, the p rotein is a 
recombinant protein; and the solution is an aqueous solution. In certain embodiments, the 
protein is a recombinant human antibody; and the solution is an aqueous solution. In certain 
embodiments, the protein is a recombinant human protein; and the solution is an aqueous 
solution. 

In a third preferred embodiment, the present invention relates to a method of decreasing 
the toxicological risk associated with administering a protein t o a mammal i n need thereof, 
comprising the steps of adding to a first solution of a protein a compound of the present 
invention to give a second solution; and administering to a mammal in need thereof a 
therapeutic amount of said second solution. In certain embodiments, the protein is a 
recombinant protein. In certain embodiments, the protein is a recombinant antibody. In 
certain embodiments, the protein is a recombinant human antibody. In certain embodiments, 
the protein is a recombinant mammalian protein. In certain embodiments, the protein is a 
recombinant human protein. In certain embodiments, the protein is recombinant human 
insulin, recombinant human erythropoietin or a recombinant human interferon. In certain 
embodiments, the first solution and the second solution are aqueous solutions. In certain 
embodiments, the protein is a recombinant protein; and the first solution and the second 
solution are aqueous solutions. In certain embodiments, the protein is a recombinant human 
antibody; and the first solution and the second solution are aqueous solutions. In certain 
embodiments, the protein is a recombinant human protein; and the first solution and the second 
solution are aqueous solutions. 

In another preferred embodiment, the present invention relates to a method of 
facilitating native folding of a recombinant protein in solution, comprising the step of 
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combining in a solution a compound of the present invention and a recombinant protein. In 
certain embodiments, the recombinant protein is a recombinant antibody. In certain 
embodiments, the recombinant protein is a recombinant human antibody. In certain 
embodiments, the recombinant protein is a recombinant mammalian protein. In certain 
embodiments, the recombinant protein is a recombinant human protein. In certain 
embodiments, the recombinant protein is recombinant human insulin, recombinant human 
erythropoietin or a recombinant human interferon. In certain embodiments, the solution is an 
aqueous solution. In certain embodiments, the recombinant protein is a recombinant human 
antibody; and the solution is an aqueous solution. In certain embodiments, the recombinant 
protein is a recombinant human protein; and the solution is an aqueous solution. 

Kinetic model approach for stabilizing proteins towards aggregation 

To see how additives affect aggregation rate, the rate constant for aggregation, kagg, can 
be expressed using transition state theory as: 



(5) 



where is Boltzmann's constant, T is the absolute temperature, h, is Planck's constant, and 
is the equilibrium constant between the reactants and the transition state for the reaction (either 
equation 2 or 4). The change in relative reaction rate due to an additive {X) at constant 

temperature and pressure can therefore be expressed as: 



dm. 



J 



'T,P,mp V X yT,P,mp 

where ntx is the molality of additive. Using the Wyman linkage relation, the above expression 
can be written in terms of the extent of binding of the additive to the protein species; 
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where Zx is the thermodynamic activity of additive, and each F is a preferential binding 
coefficient. Wyman Jr., J. Adv. Protein Chem. 1964, 19, 223-286; Timasheff, S. N. PNAS 2002, 
99, 9721-9726; Baynes, B. M.; Trout, B. L. J. Phys. Chem. B 2003, submitted for publication. 
r*px is the number of additive molecules bound to the transition state of equation 2 or 4, and 
r^px is the number of additive molecules bound to the reactant in the same equation. Since 
(51n ax I dmx)T.p.mp is positive, equation 8 shows that in order for an additive to decrease the 
rate of aggregation, the additive must bind less to the transition state than to the reactant, 
making T^xp - ^^xp negative. 

Attenuation of protein aggregation 

In the pharmaceutical industry today, a refolding buffer additive used to increase the 
yield of active protein is the amino acid L-arginine. Arginine can increase the yield of active 
protein by an order of magnitude or more and has been successfully used to refold proteins 
such as tPA, interferon y, lysozyme, carbonic anhydrase B, factor XIII, and antibodies. 
Rudolph, R.; Fischer, S.; Mattes, R. "German Patent DE 3 537 708, Process for the Activation of 
tPA after Expression in Prokaryotic Cells", 1985; Arora, D.; Khanna, N. J. Biotechnol. 1996, 52, 
127-133; Armstrong, N.; de Lencastre, A.; Gouaux, E. Protein Sci. 1999, 8, 1475-1483; Rinas, 
U.; Risse, B.; Jaenicke, R.; Abel, K.-J., Zettleneissl, G. Biol. Chem. Hoppe-Seyler 1990, 371, 49- 
56; Buchner, J. ; Rudolph, R. Biotechnology 1 991, 9, 1 57-162. A rginine has been shown to 
increase the yield of renatured protein by decreasing the rate of aggregation. Hevehan, D. L.; 
Clark, E. D. B. Biotechnol. Bioeng. 1997, 54, 221-230. It is presently disclosed that arginine 
has a critical combination of two simple factors that enable it to prevent aggregation during 
folding. These factors include size and binding. 

1 . Size. Arginine is a much larger molecule than water, the primary solvent. 

2. Binding. Protein molecules in isolation do not have a significant preference to 
be solvated by either arginine or water. 

The mechanism by which the factors above affect aggregation is shown schematically 
in Figure 1 . As the protein molecules diffuse toward each other, the size property ensures that 
a region of preferential hydration will form between the protein molecules because water but 
not the additive can fit in the gap (the oval in the transition state A2* of Figure 1). This is 
analogous to "osmotic stress" effects on the equilibrium between two macroniolecular 
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conformations where one conformation has a crevice that water can enter but an additive 
cannot. Parsegian, V. A.; Rand, R. P.; Rav, D. C. PNAS USA 2000, 97, 3987-3992. The 
binding property ensures that when there is no steric constraint due to such a gap, arginine and 
water can solvate the protein equally well. This means that the region of preferential 
hydration shown in Figure 1 is the only contribution to the preferential binding coefficients of 
the additive with the protein in any of the three states shown (U + U, A2*, A2). Because the 
transition state is preferentially hydrated, T^xp is negative. Therefore the quantity r^;o> -rVpis 
negative and aggregation is slowed. Any additive that has these two properties will deter 
aggregation during folding or in any other situation where a bimolecular step is rate limiting. 

The size and binding properties are both necessary for prevention of aggregation. 
Molecules that meet the size criterion but not the binding criterion will either accelerate 
aggregation (such as "crowders" like dextran) or be denaturants (such as guanidinium 
chloride) and therefore have other undesirable effects on protein stability. Linder, R.; Ralston, 
G. Biophys. Chem. 1995, 57; 15-25; Orsini, G.; Goldberg, M. E. J. Biol. Chem. 1978, 253, 3453- 
3458; Jasuja, R. "B-114R Protein Drugs: Manufacturing Technologies", Technical Report, 
Business Communications Company, Inc., 2000. A molecule that does not meet the size 
criterion but meets the binding criterion will have almost no effect on aggregation. 

The two properties above differentiate molecules that may have advantageous effects 
on aggregation via the mechanism above from those that may not. It is believed that there are 
many molecules that have not been used as additives which have both of the above properties. 
Since these properties are presently disclosed, arginine was not selected with them in mind, 
implying that another yet untested molecule may exemplify the properties to a larger extent 
and have superior aggregation preventing characteristics. As non-limiting examples, some 
molecules with the two properties above that may prevent aggregation via a similar 
mechanism include: 

• Citrulline 

• Arginine or citrulline derivatives with a longer or shorter methylene linker between 
the amino acid backbone and guanidino or urea group (Figure 2). 

• Arginine or citrulline derivatives where the amino acid backbone group is replaced 
by another large functional group which does not bind to proteins. (For example, 2- 
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guanidino acetic acid, 3-guanidino propanoic acid, 4-guanidino butyric acid, 5- 
guanidino pentanoic acid, etc.) 

• Molecules that are not randomly orientated in solution near proteins. Such 
molecules can be constructed by covalently attaching a molecule which stabilizes 
proteins against unfolding with a molecule that destabilizes proteins against 
unfolding. Examples of novel molecules designed based on this idea are shown in 
Figure 3. A partial list of molecules that are known to stabilize and destabilize 
proteins against unfolding are shown in Table 1. 



Table 1. 



Protein Stabilizer 


Protein Destabilizer 


Sugars (e.g. glucose, sucrose, trehalose) 


Urea 


Polyols (e.g. sorbitol, mannitol) 


Guanidinium chloride 


Dextran 


Detergents (e.g. sodium dodecyl sulfate, Tris) 


Kosmotropes 


Chaotropes 


Glycine, glycine betaine 





Compounds of the present invention 

Based on the studies described in the previous section, compounds of the present 
invention may be prepared by functionalizing a molecule that does not bind to a protein with at 
least one protein binding group. In other words, compounds of the present invention possess a 
non protein bonding moiety and a protein binding group. Molecules that do not bind to proteins 
include but are not limited to osmolytes and kosmotropes, such as glycerol, glycine betaine, 
dendrimers, and trimethyl amine N-oxide. Other such molecules are known to those skilled in 
the art. 

A protein-binding group is a molecule or functional group that binds to some proteins. 
Many molecules that fall in this class are, for example, denaturants or surfactants. Some non- 
limiting examples of protein-binding molecules are: the guanidinium ion, urea, amino acids 
(such as arginine, lysine, aspartate, glutamate), sodium dodecyl sulfate, tweens (polysorbate), 
poloxamers, and ions (such as citrate and acetate). A group or molecule does not need to bind to 
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all proteins to be classified as a "protein-binding group;" rather, it merely needs to bind to some 
proteins. The concepts of "binding" and groups or molecules that bind to proteins are well- 
known to those skilled in the art. 

The net effect of functionalizing a non-binder with a protein-binding group will be to 
move the protein preferential binding coefficient toward zero. Molecules that are large, but have 
a protein preferential binding coefficient near zero, have the properties that they prevent 
aggregation but do not destabilize native protein molecules. Thus, these molecules are useful as 

anti-aggregation additives. 

Statistical model approach for stabilizins proteins towards a£sreeation 

Additives perturb the chemical potential of the protein system by associating either more 
strongly or more weakly with the protein than water. This phenomenon, called "preferential 
binding," is of great interest because it governs the physical and chemical properties of proteins. 
Timasheff; S. N.Adv. Protein Chem. 1998, 51, 355-431. 

When an additive (X) is added to an aqueous protein solution, it alters the chemical 
potential of the protein {\ip) via the following relationship: 



where A}x/> is the transfer free energy of the protein from pure water into the mixed solvent 
system, m is molality, and subscripts A' and P identify the additive and protein respectively. Lee, 
J. C; Timasheff, S. N. J. Biol. Chem. 1981, 256, 7193-7201. Two partial derivatives appear in 
equation 10. The first captures the dependence of the additive chemical potential on additive 
molality and can be evaluated by experiments on a binary mixture of additive and water {mp — > 
0). The second partial derivative is the "preferential binding coefficient;" Txp: 




(9), (10) 
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The preferential binding coefficient is a way in which binding can be defined 
thermodynamically. It is also particularly usefUl when binding is weak. The preferential binding 
coefficient is a measure of the excess number of additive molecules in the domain of the protein 
per protein molecule (Figure 4). The connection between the thermodynamic definition 
(equation 11) and the intuitive notion of binding (local excess number of molecules) comes from 
statistical mechanics; where it can be shown that: 



In the above equation, n denotes the number of a specific type of molecule (subscript X 
for the additive and subscript W for water) in a certain domain (superscript / for a bulk volume 
outside of the vicinity of the protein and superscript // for a volume in the protein vicinity), and 
angle brackets denote an ensemble average. Kirkwood, J. G.; Goldberg, R. J. J. Chem. Plays. 
1950, 18, 54-57; Schellman, J. A. Biopolymers 1978, 17, 1305-1322. Note that Fxp is 
independent of the choice of the boundary between the domains, as long as the boundary is far 
enough from the protein. 

If the additive concentration is higher in the vicinity of the protein than in the bulk, Fxp is 
greater than zero, and \xp is lower in the presence of the additive than in its absence. Denaturants 
such as urea and guanidinium chloride exhibit this type of binding behavior. The reverse is true 
for sugars, such as trehalose. In trehalose solutions, there is generally a deficiency of trehalose 
and an excess of water in the vicinity of the protein. For this "preferential hydration" case, Fxp is 
less than zero, and |a/? is higher in the presence of the additive. 

Timasheff pioneered the use of high-precision densitometry to measure preferential 
binding coefficients for protein-cosolvent systems. Lee, J. C; Timasheff, S. N. J. Biol. Chem. 
1981, 256, 7193-7201; Lee; I. C.; Timasheff; S. N. Biochemistry 1974, 13. 257-265; Gekko, 
K.; Timasheff, S. N. Biochemistry 1981, 20. 4667-4676; Gekko, K.; Timasheff, S. N. 
Biochemistry 1981, 20, 4677-4686. More recently, differential scanning calorimetry (DSC) 
and vapor pressure osmometry (VPO) have been used to the same end. Poklar, N.; Petrovcic. N.; 
Oblak, M.; Vesnaver; G. Protein Sci. 1999, 8, 832-840; Courtenay, E. S.: Capp, M. W.; 
Anderson; C. F.; Record Jr., 11. 7. B iochemistry 2 000^ 3 9. 4455-4471. Preferential binding 
coefficients are rigorous thermodynamic quantities and are related to virial coefficients, activity 




(12) 



-25- 



coefficients, and free energies via standard thermodynamic relations for multi-component 
solutions. Casassa. E. F.; Eisenberg, Adv. Protein Chem. 1964, J 9. 287-395. 

Experimental studies by the above methods have led to some generalizations about 
preferential binding coefficients: 

1. Fxp may be positive or negative, indicating that interactions of the protein and 
additive are favorable or unfavorable, respectively. 

2. fxp is proportional to additive molality at low concentration of additive (often as high 
as /Tlx ~ 1 m and higher). Courtenay, E. S.: Capp, M. W.; Anderson; C. F.; Record 
Jr., \\. T. Biochemistry 2000, 39, 4455-4471; Greene Jr., R. F.; Pace. C. N. J. 
Biol. Chem. 1974, 249, 5388-5393; Record Jr., M. T.; Zhang; W.; Anderson; C. F. 
Adv. Protein Chem. 1998, 51, 281-353. 

3. Fxp is roughly proportional to the protein-solvent interfacial area. Lee, J. C; 
TimashefF, S. N. J. Biol. Chem. 1981, 256, 7193-7201. 

The second generalization above, together with the fact that many binary mixtures of 
additive and water (mp — > 0) are nearly ideal at low concentration of additive, leads to a useful 
simplification of equation 10: 



Equation 15 provides a simple and convenient link between preferential binding coefficients and 
free energies. This relation leads to the useful rule that when Fxp is proportional to mx, for each 
additive molecule that preferentially interacts with the protein, the protein's free energy is 
reduced by approximately 0.6 kcal/mol at 25''C. The simplicity of this relation is a natural result 
of the close relationship between Txp and a second virial coefficient. 

To be able to predict preferential binding coefficients and understand their origins, the 
above thermodynamic framework and general observations must be augmented by a mechanistic 
model. Several such models have been presented in the literature, including models based on the 
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binding polynomial or statistical mechanical partition function, solvent-additive exchange at 
defined sites, additive partitioning between the local and bulk domains, and group contribution 
methods for estimating transfer fi-ee energies. 

The most general model of additive binding hitherto presented comes from considering 
an equilibrium of all possible protein-additive complexes, fi-om which it can be shown that: 

/^^i; =-RT\n{\ + YL^ym'„mi) (16) 
'■ J 

where Ky is the equilibrium constant for a reaction of a protein molecule, i molecules of water, 
and j molecules of additive into a complex. Wyman, J.; Gill; S. J. Binding and Linkage: 
Functional Chemistry of Biological Macromolecules: University Science Books: 1990. 
While this model is completely general, its utility is limited because it is not possible to 
determine experimentally the many Ky parameters present in equation 16. 

Schellman's site exchange model, provides a way to simplify this general expression to a 
form containing a single parameter. Schellman, J. A. Biopolymers 1978, 77, 1305-1322. This 
model treats binding as a family of protein-solvent exchange reactions such as: 

PWi + X^PX + iW (17) 

where P is the protein, W is water, X is cosolvent; and i is the exchange stoichiometry. The 
simplification requires the assumptions that 1:1 exchange reactions (i = 1) occur on a fixed 
number of identical, independent sites and that the sites are far fi"om saturation with additive (i.e. 
the apparent dissociation equilibrium constant for each site is well above the additive 
concentration). The number of sites, n, is approximated by the number of water molecules 
present in a monolayer around the protein. These simplifications reduce equation 16 to: 

A^'; =-nRT{K)mj, (18) 

where {fO is the average equilibrium constant of binding at a single site. The single parameter 
{K) can then be determined fi"om an experimental measurement of Fxp. When equation 15 holds, 
the relation between {K) and Fxp is simply: 

{K) = T^,lnm, (19) 
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Values of (lO for different proteins in this linear regime are roughly equal. Schellman, J. A. 
Biophys. Chem. 2002, 96. 91-101. (/O cannot, however, be determined without knowledge of 
Fxp or other free energy data on the particular additive system of interest. In fact, one can say 

that (/O is defined by Fxp. 

Another model that recasts preferential binding coefficient data in terms of a single 
model parameter is the local-bulk domain model developed by Courtenay et al . Courtenay, E. 
S.: Capp, M. W.; Anderson; C. F.; Record Jr., 11. T. Biochemistry 2000, 39, 4455-4471. The 
parameter in this model is the partition coefficient K^, relating the number of water molecules 
and additive molecules in the local and bulk domains via: 



Similar to the site exchange model, the convention used in this model is that the local domain 
consists of a monolayer of water and enough additive to obtain the experimentally observed Fxp. 
Note that because the absolute occupancy of water and additive in the local domain carmot be 
easily determined by experiment, the local-bulk domain model effectively defines nw. Like {fO, 
values of Kp can be used to predict Fxp at other additive concentrations or for other proteins in 
the same additive, but predictions cannot be made in the absence of Fxp or free energy data on 
the same additive system. 

Lastly, transfer free energy models, pioneered by Bolen's group, take a different 
approach. Liu, Y.F.; Bolen, D. ^.Biochemistry 1995, 34. 12884-12891. T hese models 
conceptually divide whole proteins into groups such as the amino acid side chains and the 
protein backbone and model the transfer free energy of the whole protein as a sum of the transfer 
free energy of the groups it comprises, via: 



where Ag, is the transfer free energy of the model group and a\ is the solvent accessible area of 
the group in the whole protein, normalized to the solvent accessible area of the model compound. 
Tanford, C. J. Am. Chem. Soc. 1964, 86, 2050-2059. The overall Afx"> can then be predicted 
for any system of known structure. In the context of the previously described models, the 
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transfer free energy model can be thought of as a Hnearized binding model where each surface 
group or amino acid in the protein represents a different type of independent binding site, and the 
binding constants for those sites are determined by experiments on model compounds, such as 
free amino acids or cyclic di-amino acid compounds. Predictions made by transfer free energy 
models have met with mixed success. A linear group contribution model (equation 21) may be 
too simple to capture all of the important contributions to A^"p• Bolen, D. W. Protein 
Stabilizaiton by Naturally Occurring Osmolytes. In Protein Structure, Stability, and 
Folding; Humana Press: 2001. 

While the above models have helped in the understanding of the phenomenon of 
preferential binding, they generally incorporate strong assumptions, and they necessitate the use 
of experimental data on highly analogous systems in order to determine model parameters and 
make predictions. Thus, their uses as predictive tools and as tools to gain insight into specific 
systems are limited. 

One aspect of the present invention relates to a predictive, molecular-level approach for 
the s tudy o f p referential b inding b ased o n all-atom, s tatistical m echanical m odels t hat u se n o 
adjustable parameters. To date, statistical mechanical models of preferential binding have only 
been developed for interactions of ions with charged cylinders and for interactions of two- 
dimensional, "hard circles" with a linear interface, both far too simple to be generally applied to 
protein-additive systems. Anderson; C. F.; Record Jr., M. T. J. Phys. Chem. 1993, 97, 7116- 
7126; Mills, P.; Anderson, C. F.; Record Jr., M. T. J. Phys. Chem. 1986, PO, 6 541-6548; 
Tang. K. E. S.: Bloomfield, V. A. Biophys. J. 2002, 82. 2876-2991. Other explicit mixed 
solvent simulations of proteins and amino acids have been performed, but these studies did not 
compute thermodynamic quantities related to preferential binding. Zou, Q.; Bennion. B . J.; 
Daggett, v.; Murphy, K. P. J. Am. Chem. Soc. 2002, 124, 1192-1202; Bennion, B. J.; 
Daggett, V . PA'/^^ 2 003, 100. 5142-5147; Tirado-Rives, J.;Orozco, M.; Jorgensen, W. L. 
Biochemistry 1997, 36, 7313-7329; Alonso, D. O. V.; Daggett, V. J. Mol. Biol. 1995, 247, 
501-520; Caflisch. A.;Karplus, XI. Structt. Fold. Des. 1999, 7. 477-488. In the present 
invention, the number of "bound" molecules are defined in a thermodynamically consistent way 
and do not a priori incorporate any information about "binding sites." The use of this approach 
for the computation of preferential binding coefficients was validated in two systems by 
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comparison with experimental data from the Hterature. Additionally, the molecular-level detail 
of the approach provides new insights into the following issues: 

1 . The changes in solvent and additive concentration as a function of distance from the 
protein surface. 

2. A precise definition of the "local domain" (Figure 4). 

3. The differences in preferential binding or apparent binding equilibrium constant at 
different locations on the protein-solvent interface. 

The success of this method in modeling preferential binding indicates that it captures the 
important underlying physics of protein-additive-water systems and that the difficulty in 
quantitative prediction to date can be surmounted by explicitly incorporating the complex 
protein-solvent and solvent-solvent interactions. 

A Molecular-Level Approach to Computins Preferential Binding 

One aspect of the present invention relates to the use of explicit atomic interaction 
potentials (force fields), such as Lennard- Jones, Coulombic, spring, and torsion interactions, with 
pre-fit coefficients. Brooks; B. R.; Bruccoleri; R. E.; Olafson, B. D.; States, D. J.; 
Swaminathan, W.: Karplus, M. J. Comp. Chem. 1983, 4, 187-217; Ha; S. N.; Giammona; 
A.: Field, M.; Brady, J. W. Carbohydrate Res. 1988, 180, 207-221. Thermodynamic 
properties, such as preferential binding coefficients, are computed by averaging in the time 
domain via molecular dynamics (MD). A snapshot from a dynamic simulation of RNase Tl in a 
urea solution is shown in Figure 5, which was generated with VMD. Humphrey, W.; Dalke, A.; 
Schulten, K. J. Molec. Graphics 1996, 14, 33-38. The results of the simulations contain all of 
the information needed to extract thermodynamic properties, such as Fxp. 

Molecular dynamics uses Newton's second law of motion, that acceleration is the 
quotient of force and mass, to compute the positions of each atom in the system as a frinction of 
time. To do this, an energy model, sometimes called a "force field," that can be used to compute 
the net force on any atom in any configuration is employed. 

During the MD run, the positions of each atom are recorded at fixed intervals in time. 
These "snapshots" form an ensemble of configurations which can then be used to compute 
thermodynamic properties, such as Fxp. 
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Importantly, this method of computing Fxp does not introduce any adjustable parameters 
to model preferential binding or any other aspect of a system containing a protein and solvent- 
additive components. All of parameters required by the MD method for energy computations are 
determined independently of this particular modeling objective, and in fact have been shown to 
be generally applicable to biological systems. Karplus, M., McCammon, J. A. Nature. Struct. 
Biol. 2002, 9, 646-652. Thus, the method developed here could be used to estimate Fxp and 
A^")? in systems where no experimental data is available. It therefore facilitates the study of 
preferential binding when direct experimental study is difficult, such as at transition state 
configurations or at marginally stable states of proteins. Furthermore, it yields detailed, local, 
molecular-level insight into the system studied. 

Another benefit of this approach is that when equation 15 holds (such as for urea and 
glycerol), the protein transfer fi-ee energy (Aji"p) can be calculated from a single Fxp simulation. 
Traditional free energy calculation methods such as thermodynamic integration require 15-20 

trajectories, which is computationally difficult for protein systems of this size. Bash, P. A.; 
Singh, U. C: Langridge, R..; Kollman. P. A. Science 87, 236. 564-569; Kollman, P. Chem. 
Rev. 1993, P3. 2395-2417. 

Preferential Bindins Coefficients of Constituent Groups 

Because proteins have a range of different functional groups in different orientations on 
their surfaces, the concentrations of solvents and additives near different patches on the protein's 
surface may be different. For example, the vicinity of a hydrophobic patch on the protein may 
have a lower concentration of water and a higher concentration of additive than in the vicinity of 
a hydrophilic patch. Preferential binding experiments capture only the average effect arising 
from all of the interactions over the entire protein-solvent interface; however, molecular 
simulations allow more detailed analyses of the local contributions to preferential binding 
coefficients. 

A protein can be thought of as a set of non-overlapping constituent groups, each of which 
has its own preferential binding coefficient defined by the composition of the solvent in its 
immediate vicinity. Tanford, C. J. Am. Chem. Soc. 1964, 86, 2050-2059. Similar to group 
contribution methods for computing transfer fi-ee energies, one possible group definition is that 
each type of amino acid side chain (up to 20) and the amino acid backbone are distinct groups. 
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To compute a preferential binding coefficient for a constituent group, the solvent molecules in 
the local domain are assigned only to the nearest group (i), and the "group preferential binding 
coefficients" (Fxp, i) can be defined as: 



where n x,/ and n w,/ are the number of additive and water molecules in the local domain that are 
nearest to group /. If each additive molecule in the local domain is assigned to a group, the 
overall preferential binding coefficient is simply the sum of all of the group preferential binding 
coefficients: 



The group preferential binding coefficients decompose the effect of each small subset of 
the protein on the overall preferential binding coefficient. This is analogous to the group 
contribution models for transfer fi-ee energy except that the parameters are extracted from a 
simulation of an entire protein instead of experiments on model compounds. 

Minimum Simulation Time 

Sufficient sampling of position-space configurations in time is required for the accurate 
calculation of Fxp via equation 1 1 . Assuming that the average protein solution structure is close 
to that of the initial (crystal) structure and that water molecules sample position space rapidly 
because of their high density, the most important time scale to be captured is that of the additives 
sampling position space. One way to estimate this time is that it must be much larger than the 
average time between additive-additive contacts. 

An estimate of the time between contacts can be obtained as: 
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where D is the additive diffiisivity, Wsoiv is the solvent volume, and nx is the number of additive 
molecules. For the simulations performed here, the solvent is mostly water, so equation 24 can 
be further simplified to yield: 
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where Na is Avogadro's number and pw is the density of water in kg/m''. For aim additive in 
water system with a additive diffusivity of 2x10'' mVs (a lower bound on the diffusivities of the 
additives studied here), Xcomaa is about 30 ps. Thus, nanosecond trajectories will be required for 
good sampling of additive position space. Importantly, this time increases as the additive 
concentration decreases, implying that there is a minimum concentration that can be studied with 
any given amount of computational resources. 

Radial Distribution Functions of Water and Additives 

The radial distribution functions of water, urea, and glycerol were computed for all three 
simulations as described in the Exemplification section and are shown in Figure 6. 

At very short distances, r < 0.6A for water and r < l.OA for glycerol and urea, regions of 
total solvent and additive exclusion due to very strong van der Waals repulsion can be seen. The 
size oft hese "totally excluded" r egions i s m uch s mailer t han o ne w ould e xpect b ased o n t he 
apparent van der Waals radii of the solvent and additive molecules alone (for example, r w l.SA 
for water and 2.2A for urea), indicating that electrostatic attractive forces play an important role 
in solvation even at these distances. Schellman, J. A. Biophys. J. 2003, 85. 108-125. After 
the regions of total exclusion, strong first coordination shells of these three molecules can be 
clearly seen. The peaks of the first coordination shells become more distant fi-om the protein as 
the size of the molecules they correspond to increases. Significantly smaller second coordination 
shell peaks are also visible for urea solvating RNase Tl and glycerol solvating RNase A. At 
distances greater than 6-7A fi'om the protein, solvation shells cannot be discerned, and the 
number densities of water, urea, and glycerol reach their bulk values. 

In the simulations of RNase Tl in glycerol and urea solutions, the radial distribution 
functions for glycerol and urea are quite different. The maximum value of gx(r) for urea is over 
4.5, while that for glycerol is about 2.5. The difference in these maximum values, while 
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significant, is not sufficient to say that the number of urea molecules coordinated to the protein 
(rix) is higher than the number of glycerol molecules coordinated, this can only be done by 
integrating each gx{r) function appropriately via equation 31. 

The radial distribution functions for both water and glycerol are similar in the simulations 
of RNase A and RNase Tl in glycerol solution, despite the fact that the proteins and the pHs of 
the solutions are different. Given that the proteins are of similar size, this observation is 
consistent with the fact that the values of Fxp for the two solutions are close. 

Preferential Binding Coefficients 

The radial distribution functions in Figure 6 suggest that r* in the range of 6-8A is an 
appropriate choice of boundary between the local and bulk domains. The error in Fxp introduced 
by a particular choice of the boundary distance, r», can be estimated by plotting the apparent 
preferential binding coefficient (Fxp) versus r* (Figure 7). Fxp depends very strongly on r» in the 
first solvation shell (r = 0 - 4A) and weakly on r* in the second solvation shell (/• = 4 - 6A). In 
the range r = 6 - SA, the dependence of Fxp on r* is small (±0.5), and is less than the statistical 
error in Fxp (shown in Table 2, explained below). Therefore, a cutoff distance of 6A, or about 
two solvation shells, is sufficiently large to minimize systematic error in Fxp caused by the choice 
of r*. If only a single solvation shell were considered (r* ~ 3.5 - 4A), a systematic error in Fxp of 
approximately 0.5 - 1 molecules would be introduced as a result of neglect of the second 
solvation shell. 

The preferential binding coefficient, Fxp, was computed via equation 1 1 using r* = 6A as 
the boundary between the local and bulk domains. A confidence interval for this ensemble 
average was computed as described in the Exemplification section. The binding coefficients and 
their statistical uncertainties are shown in Table 2. 



Table 2. Preferential binding coefficients computed fi-om MD simulations and compared with 
available experimental data at similar additive concentrations. 



System 


mbiiik 


Simulation Fxp 


Experimental Fxp 


Urea / Rnase Tl 


1.10m 


5.2 ± 1.0 


6.4^ 


Glycerol / Rnase Tl 


1.07 m 


-1.6 ±0.8 
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Glycerol / Rnase A 


0.91 m 


-0.9 ± 1.0 


-1.7 ±0.8" 


" Lin, T. Y.; Timasheff, S. N. Biochemistry 94 


, 33, 12695-12701. 



Gekko, K.; Timasheff, S. N. Biochemistry 1981, 20. 4667-4676. 



A wide range of behavior (positive and negative preferential binding coefficients) can be 
modeled without the use of adjustable parameters. The confidence intervals on FxpCMD) are an 
estimate of the statistical error resulting from the use of a finite trajectory. For easier 

comparison, the experimental values of Fxp reported above were interpolated to mt,uik from data 
sets spanning the molality of interest. 

Experimental values from the literature were available for two out of three of these 
protein-additive systems, and the computed values of Fxp agree quite favorably with these values. 
The fact that this occurs for both positive and negative values of Fxp without the use of any 
adjustable parameters is very encouraging. For an additive that obeys equation 15, the 
confidence intervals of ±1.0 in Fxp represents a confidence limit in the transfer free energy of 
about 0.6 kcal/mol, which is a typical value for free energies calculated via this type of molecular 
simulation. Achievement of this level of accuracy despite the fact that structural fluctuations in 
the native state ensemble of proteins have been observed on much longer time scales than the 
time scale of the simulations performed here suggests that solvent dynamics are more important 
than protein structural dynamics in determining Fxp. Duan, Y.; Kollman, P. A. Science 1998, 
282, 740-744. 

Fxp(0 probability density functions for the simulations of RNase Tl in urea and glycerol 
solution are shown in Figure 8. The range of instantaneous values of the preferential binding 
coefficient, Fxp(0, is quite large relative to the absolute values of Fxp. FxpCO values in excess of 
Fxp ± 15 are observed. The breadths of these distributions are related to the size of the interface 
between the local and bulk domains and indicate the importance of sampling a large number of 
solvent configurations to obtain the macroscopic, averaged Fxp (equation 27). 

TTie Relation between Solvent Accessible Area and the Number of Molecules in the Local 
Domain 

The solvent accessible areas of whole proteins (S AA) and constituent groups (SAA,) in 
crystal structures have been used extensively in analyzing proteins. SAA and SAA, are 
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essentially simple ways of measuring water coordination numbers. In models developed to date, 
SAA or SAA„ has been used to estimate /i,v or Mw,; by assuming that the local domain is a 
inonolayer of water and each water molecule occupies approximately lOA^ of the solvent 

accessible area. Since the present invention introduces a new notion of the local domain, it is 
worthwhile to see what relationships exist between SAAj and the coordination numbers «w,/ and 
nx,i that utilize this defmition. 

A scatter plot of the solvent accessible area of a set of constituent groups (amino acid side 
chains and the protein backbone) versus the number of water molecules in the local domain for 
three different simulations is shown in Figure 9. Solvent accessible area was calculated 
analytically in CHARMM (based on Richmond's method) using a 1 .4A probe. Richmond, T. J. 
J. Mol. Biol. 1 984, 178, 63-89. There is a strong, linear correlation of these variables with 
slope 4.2 AVmolecule and correlation coefficient 0.96. Similarly strong correlations are seen for 
SAA, with /Zx,! in individual simulations. A summary of proportionality constants and correlation 
coefficients for these relationships is shown in Table 3. If the time average SAA, from each 
dynamics simulation is used instead of the crystal structure SAA, values, the correlation 
coefficients increase slightly. Because the time average solvent accessible areas are higher than 
those in the crystal structure, the proportionality constants shown in Table 3 also increase. 



Table 3. Relationships between solvent accessible area in each protein crystal structure and 
number of solvent molecules in the local domain for different protein-additive systems, 
symbolizes the correlation coefficient. 



Species (/) 


Protein 


Avg. Protein SAA/n/' 
(AVmolecule) 




Water 


RNase A/Tl 


4.2 


0.96 


0.91 m Glycerol 


RNase A 


290 


0.96 


1.07 m Glycerol 


RNase Tl 


230 


0.93 


1.10 m Glycerol 


RNase Tl 


170 


0.98 



Constituent Group Preferential Binding Coefficients 

The constituent group preferential binding coefficients were calculated for each 
simulation as described in the Exemplification section and are shown in Figures 10 - 13 as the 
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number of water and additive molecules coordinated to each constituent group. In each figure, a 
line at the bulk solution composition is also plotted, enabling a quick determination of the 
composition of the solvent in the vicinity of a constituent group compared to the bulk solvent. 
The statistical uncertainties in the values of and n'^^j (and consequently Fxp,,) are high. 
Because of these uncertainties, we will not report specific values of the group preferential 
binding coefficients, but rather classify them into broad categories based on their statistical 
likelihood of being either positive, negative, or zero/ indeterminate. 

The average number of water and glycerol molecules coordinated to each of the 15 serine 
residues in RNase Tl are shown in Figure 10. A wide range of binding behavior can be seen 
among the serine residues, all of which have a good degree of solvent exposure. Ser 17, 35, and 
72 fall above the bulk concentration line and have positive preferential binding coefficients, Ser 
63 falls below the line and has a negative preferential binding coefficient, and the preferential 
binding coefficients of the remaining 1 1 serine residues are not statistically different from zero. 
The wide range of local concentrations in the vicinities of these serine residues indicates that 
developing a group contribution method to estimate Fxp or A/u'p based on primary sequence 
information and solvent accessibility (n^'w,/) alone may be difficult. In addition to the type of 
amino acids present at the protein-solvent interface, other effects such as specific combinations 
of residues and secondary or tertiary structure must be important in determining water and 
additive binding behavior. These factors probably contribute to the range of local concentrations 
seen in Figure 10. For example, Ser35 and Ser72 are proximal to each other and several Gly and 
Tyr side chains (Gly 34, 70, 71, and Tyr 68), which tend to have positive preferential binding 
coefficients in- glycerol (Figure 12). This may be the reason that the group preferential binding 
coefficients for these residues are higher than those of the other serine residues. 

The preferential binding behavior of urea and glycerol, with each type of amino acid in 
RNase Tl and the protein backbone are shown in Figures 1 1 and 12. In urea solution, the protein 
backbone and Ser as well as the hydrophobic amino acid side chains of Cys, Gly, Len, Phe, Pro, 
Tyr, and Val all preferentially bind urea, while the hydrophilic Asp preferentially binds water. In 
glycerol solution, only Tyr and Gly preferentially bind glycerol, and Asp and Glu preferentially 
bind water. Qualitatively, the binding behavior of the amino acid side chains of RNase Tl 
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follow a hydrophobic series, with the hydrophobic side chains tending to bind more additive and 
the hydrophilic ones tending to bind more water. 

The binding behavior of glycerol and water with the amino acid side chains and backbone 
in RNase A, shown in Figure 13, is significantly d ifferent than the binding behavior of these 
solvent components with the same constituent groups in RNase Tl. (Note that the protonation 
states o f A sp, G lu, a nd His a re d ifferent i n t he two s imulations.) T he a mino a cid b ackbone, 
which occupies a large fraction of the protein-solvent interface as indicated by its high value of 
n"^^i, has a binding coefficient near zero in RNase Tl and a significant negative binding 
coefficient in RNase A. More strikingly, Tyr in RNase Tl preferentially binds glycerol whereas 
Tyr in RNase A preferentially binds water. This is likely because the six Tyr residues in RNase 
A are at or near the solvent interface (a more hydrophilic region) whereas the nine in RNase Tl 
are mostly buried (a more hydrophobic region). This difference in solvent exposure is evident 
from the crystal structures of the proteins but also can be discerned by comparing the water 
coordination numbers for Tyr in the two proteins: n"v,j for Tyr in RNase A is higher than in 
RNase Tl, even though there are 50% more Tyr residues in RNase Tl. 

Based on the above observations, some generalizations about the effects that these 
additives have on protein folding equilibria can be postulated, the validity of which must be 
confirmed via fijture studies. In urea solution, most of the constituent groups in RNase Tl either 
preferentially bind urea or are indifferent to urea and water. Asp, which is found on the surface 
of RNase Tl, is the only constituent group that is significantly below the bulk concentration line 
in Figure 11 and therefore preferentially binds water over urea. Since the amino acids that 
compose the core of RNase Tl and are exposed upon unfolding preferentially bind urea, this 
pattern suggests that the preferential binding coefficient or urea with unfolded RNase Tl is 
higher than that with native RNase Tl. This is thermodynamically consistent with urea's well- 
known ability as a denaturant. Inversely, in glycerol solution, almost all of the constituent 
groups in RNase A and Tl are neutral or preferentially bind water. This is consistent with the 
fact that glycerol binds less to the unfolded protein than the native state, and therefore is a 
protein stabilizer. Both of these generalizations are consistent with earlier work on model 
compounds. Bolen, D. W. Protein Stabilizaiton by Naturally Occurring Osmolytes. In 
Protein Structure, Stability, and Folding; Humana Press: 2001. 
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Exemplification 

The invention now being generally described, it will be more readily understood by 
reference to the following examples, which are included merely for purposes of illustration of 
certain aspects and embodiments of the present invention, and are not intended to limit the 
invention. 

Example 1 

Molecular Simulations - Molecular dynamics was used to sample the phase space of proteins 
solvated by water and an additive. Version 28 of the CHARMM molecular dynamics package 
was used for all simulations. Brooks; B. R.; Bruccoleri; R. E.; Olafson, B. D.; States, D. J.; 
Swaminathan, W.: Karplus, M. J. Comp. Chem. 1983, 4. 187-217. The CHARMM force- 
field was used for the protein, and the TIP3P model [32] was used for water. Jorgensen, W. L.; 
Chandrasekhar. J.; Madura, J. D.; Impey, R. W.; Klein, M. L. J. Chem. Phys. 1983, 79. 
926-935. A force-field was constructed for glycerol using the standard CHARA-IIVI 
geometries and partial charges for the atoms in a -CHOH- unit. Brooks; B. R.; Bruccoleri; R. 
E.; Olafson, B. D.; States, D. J.; Swaminathan, W.: Karplus, M. J. Comp. Chem. 1983, 4, 
187-217; Ha; S. N.; Giammona; A.: Field, M.; Brady, J. W. Carbohydrate Res. 1988, ]80, 
207-221 . Urea was assumed to be planar with bond lengths equal to the CHARMM standards 
and partial charges recomputed as done previously [33] but using the CHARMM van der Waals 
mixing rules in the objective function. Duffy. E. M.; Severance. D. L., Jorgensen, W.L. Israel 
J. Chem. 1993, 33. 323-330. 

The structures of RNase A (PDB code: lfs3) and RNase Tl (PDB code: lygw) were 
obtained from the Protein Data Bank. Berman, H. M.; Westbrook, J.; Feng, Z.; Gilliand; G.; 
Bhat; T. N.; Weissig, H.; Shindyalov. I. N.; Bourne, P. E. Nucleic Acids Res. 2000, 28. 
235-242. In total; three simulations were performed: RNase A in Im glycerol (pH 3), RNase Tl 
in Im glycerol (pH 7), and RNase Tl in Im urea (pH 7). Details of each simulation are shown in 
Table 4. Each protein was solvated in a truncated octahedral box extending a minimum of 9A 
from the protein. The pH of each simulation was fixed by setting the protonation states of each 
ionizable side chain to the dominant form expected for each amino acid at the pH of interest. 
Arginine, cysteine, lysine, and tyrosine were protonated in all of the simulations. Aspartate, 
glutamate, and histidine were assumed to have pKa values of 3.4, 4.1, and 6.6, respectively; and 
were therefore protonated in the simulation at pH 3 and deprotonated at pH 7. Forsyth, W. R.; 
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Antosiewicz. J. hi.; Robertson, A. D. Proteins 2002, 48. 388-403; Edgecomb, S. P.; 
Murphy, K. P. Proteins 2002, 49, 1-6. Initial placement of water and additive molecules were 
random. Protein counterions were placed using SOLVATE 1.0. The system was first energy 
minimized at 0 K, next heated to 298.15 K, and then equilibrated for 1 nanosecond in the NTP 
ensemble at one atmosphere. For the computation of the properties of interest, two nanoseconds 
of dynamics were then run, during which statistics were computed from snapshots of the 
trajectory every picosecond. 



Table 4. Details of four molecular dynamics (AID) simulations performed, nx is the number of 
additive molecules, «w is the number of water molecules, and </> is the average dimension of the 
primary unit cell (which varies during the run at constant pressure). 



Additive 


Protein 


T(°C) 


pH 


rix 




</>(A) 


Urea 


RNaseTl 


25 


7 


90 


4274 


57.48 


Glycerol 


RNase Tl 


25 


7 


87 


4582 


59.24 


Glycerol 


RNase Tl 


25 


3 


90 


5480 


62.86 



Example 2 

Calculation of Preferential Binding Coefficients - The trajectories were then used to define the 
local and bulk regions and compute Fxp in the following manner. For the purpose of computing 
Fxp and other thermodynamic and structural parameters, each water and additive molecule was 
treated as a point at its center of mass. The distance of each of these points to the protein's van 
der Waals surface was computed, and then pvv(r) and px(A-), defined as the number densities of 
these points at a distance r from the protein, were computed. In all cases, the p(r) fimctions 
exhibited peaks and valleys characteristic of solvation shells in the range 0 < r < 6k. At 
distances in the range of 6-8A and higher, such variations are no longer seen, and the local 
number density is defined as bulk number density, p (oo). Such a region far from the protein 
containing a spatially uniform concentration of water and additive must be present in the 
simulation cell in order to define the local and bulk regions and calculate Fxp. 

The position of the boundary between the local and bulk domains, a distance of r* away 
from the surface of the protein, was then determined by choosing the minimum distance at which 
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no significant difference between p(r*) and p(oo) was apparent for either water or additive. All 
solvent molecules whose centers of mass fell inside a distance of r* from the protein's van der 
Waals surface were defined as belonging to the local domain (II), and all other solvent molecules 

were defined as belonging to the bulk domain (I). With these definitions of the domains, the 
instantaneous preferential binding coefficient, r^pit), was computed as 



for each time point in each trajectory. The preferential binding coefficient, Fxp, was then 
computed for each trajectory as the time average of these instantaneous values: 



=]]^xAndt' (27) 



The radial distribution functions gx(r) and gw(r) are defined as: 

g,(/-) = AW/A(<«) (28) 

where / represents water (ff^ or an additive (X) species. These fiinctions provide another route to 
compute Fxp: 



X 



J 



= Px (<») \sxdV - f^^T V ('"^ (29), (30), (31) 

where each integral is over the local domain or the entire system (since gx - = 0 in the bulk 
domain). 

The boundary between domains I and II must be placed far enough from the protein to 
ensure that it is in the bulk, yet at the smallest such distance so that statistical fluctuations in the 
number of molecules in the domains can be minimized. One can use the values of gjc(r) and 
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gw(r) to determine the optimal boundary. Defining Txp as the apparent preferential binding 
coefficient resulting from defining the local domain as those molecules whose centers of mass lie 
inside a distance r* fi-om the protein: 

^xpi''') = Px ksx -Sw)-rdr (32) 
The error in Fxp, Er, introduced by selecting a particular value of r* is then 



F =r* -r 

^ XP ^ XP 

dr 



. "r, ^dV , (33), (34) 



When r* is selected properly, the surface defined by r = r* is entirely in the bulk solution, gx(r*) 
= gw{r*) = 1, and Ep = 0. Thus, selecting r* as the minimum distance for which all r ^ r* satisfy 
gx{r) = gw(r) = 1 (within the error of the simulation) is optimal. 

Example 3 

Calculation of Constituent Group Preferential Binding Coefficients - For each simulation, up to 
21 constituent group preferential binding coefficients were calculated. The 21 groups were each 
type of amino acid side chain present in the protein (up to 20) and the protein backbone. The 
"protein backbone" was defined as the -NH-CH-COO- unit, as well as the two extra protons at 
the N-terminus and extra oxygen atom at the C-terminus of the protein. The glycine side chain 
was defined as the proton bound to the alpha carbon that would be replaced by a substituent to 
form a different L-amino acid. 

For the simulation of RNase Tl in glycerol solution, the constituent group preferential 
binding coefficients for the 1 5 individual serine residues in the protein were also calculated. For 
this calculation, solvent and additive molecules that were nearest to an atom in the protein that 
was not part of a serine side chain were not considered. 

Water and additive molecules were associated with a specific constituent group by 
computing the distance fi'om the center of mass of each solvent molecule to the van der Waals 
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surface of every atom in the protein, selecting the protein atom that was nearest to the solvent 
molecule, and then determining to what constituent group this nearest protein atom belonged. 

Example 4 

Estimation of Statistical Error - The statistical error arising from computing averaged properties 
from a finite trajectory was estimated in the following fashion: 

1 . The dynamic trajectory of interest was divided into n pieces. 

2. The mean of the property of interest was computed in each piece. These means were 

designated z, where i = \...n. 

3. The standard deviation of the z,- values was computed. 

4. This standard deviation was divided by n and the quotient was designated a^, an 
estimate of the error in the mean determined by time averaging the full trajectory. 

The number of pieces n into which the trajectory is divided must be small enough to ensure that 
the means of each piece (the z,) are statistically independent. An autocorrelation analysis (not 
shown) of several trajectories of rxp(0 data and the underlying molecular counts («, and «,) 
indicates that a window of about 0.2 ns is sufficiently large for this to be true. Therefore, for a 2 
ns dynamics trajectory, a value of n = 2/0.2 = 10 was used. 

For long trajectories, the statistical error a,„ is roughly proportional to the inverse square 
root of the trajectory length. This property can be used to estimate the trajectory length required 
to achieve a given level of statistical accuracy after a small trajectory has been generated and 
analyzed. 
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Incorporation by Reference 
All of the patents and publications cited herein are hereby incorporated by reference. 

Equivalents 

Those skilled in the art will recognize, or be able to ascertain using no more than routine 
experimentation, many equivalents to the specific embodiments of the invention described 
herein. Such equivalents are intended to be encompassed by the following claims. 
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We claim: 

1 . A compound, comprising a non-protein-binding moiety (NPBM) and at least one protein- 
binding group (PBG). 

2. The compound of claim 1 , wherein the NPBM is a polyol, sugar, amino acid, or 
dendrimer moiety. 

3. The compound of claim 1, wherein the PBG is a urea, guanidinium ion, detergent, amino 
acid, denaturant, surfactant, polysorbate, polaxamer, citrate, chaotrope, or acetate group. 

4. The compound of claim 2, wherein the polyol moiety is a sorbitol or mannitol moiety. 

5. The compound of claim 2, wherein the sugar moiety is a glucose, sucrose, or trehalose 
moiety. 

6. The compound of claim 2, wherein the amino acid moiety is an arginine betaine, proline, 
or ectoine moiety. 

7. The compound of claim 2, wherein the dendrimer moiety is based on benzene, 
pentaerythritol, P(CH20H)3, or TRIS. 

8. The compound of claim 3, wherein the PBG is a guanidinium ion. 

9. The compound of claim 3, wherein the PBG is sodium dodecyl sulfate. 

10. The compound of claim 1, wherein the compound has formula I: 




wherein: 

R is an electron pair, H, alkyl, aryl, heteroaryl, aralkyl, heteroaralkyl, or an alkali metal; 

R' is H, alkyl, aryl, heteroaryl, aralkyl, heteroaralkyl, or R"H2N; 

R" is an electron pair, H, alkyl, aryl, heteroaryl, aralkyl, or heteroaralkyl; 

W is O, NH2^(halogen)", or S; and 
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nis 1,2, or 4-100. 




11. 


The compound of claim 10, wherein R is an electron pair. 




12. 


The compound of claim 10, wherein R' is H. 




13. 


The compound of claim 10, wherein R' is R"H2N. 




14. 


The compound of claim 10, wherein R' is H3N*. 




15. 


The compound of claim 1 0, wherein W is NHa^Cl". 




16. 


The compound of claim 10, wherein n is 1. 




17. 


The compound of claim 10, wherein n is 2. 




18. 


The compound of claim 10, wherein n is 4. 




19. 


The compound of claim 10, wherein n is 5. 




20. 


The compound of claim 10, wherein n is 6. 




21. 


The compound of claim 10, wherein R is an electron pair, R' 


is Haisr', W is NH2^Cr, and 


n is 1. 






22. 


The compound of claim 10, wherein R is an electron pair, R' 


is HslsT, W is NHa^Cr, and 


nis2. 






23. 


The compound of claim 10, wherein R is an electron pair, R' 


is HsN^, W is NH2''C1', and 


n is 4. 






24. 


The compound of claim 10, wherein R is an electron pair, R' 


isHslsT, WisNH2^Cr,and 


n is 5. 






25. 


The compound of claim 10, wherein R is an electron pair, R' 


is HalsT, W is NH2''Cr, and 


n is 6. 






26. 


The compound of claim 10, wherein R is an electron pair, R' 


is W is 0, and nisi. 


27. 


The compound of claim 10, wherein R is an electron pair, R' 


is HaN*, W is 0, and n is 2. 


28. 


The compound of claim 10, wherein R is an electron pair, R' 


is HsN^, W is 0, and n is 4. 


29. 


The compound of claim 10, wherein R is an electron pair, R' 


is HalsT, W is 0, and n is 5. 
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30. The compound of claim 10, wherein R is an electron pair, R' is H3N*, W is O, and n is 6. 

3 1 . The compound of claim 1 0, wherein R is an electron pair, R' is H, W is NH2^Cr, and n is 
1. 

32. The compound of claim 10, wherein R is an electron pair, R' is H, W is NH2*Cr, and n is 
2. 

33. The compound of claim 10, wherein R is an electron pair, R' is H*, W is NH2 CI', and n 
is 4. 

34. The compound of claim 1 0, wherein R is an electron pair, R' is H, W is NH2"^Cr, and n is 
5. 

35. The compound of claim 10, wherein R is an electron pair, R' is H, W is NH2'^Cr, and n is 
6. 

36. The compound of claim 10, wherein R is an electron pair, R' is H, W is O, and n is 1. 

37. The compound of claim 10, wherein R is an electron pair, R' is H, W is O, and n is 2. 

38. The compound of claim 1 0, wherein R is an electron pair, R' is H, W is 0, and n is 4. 

39. The compound of claim 10, wherein R is an electron pair, R' is H, W is O, and n is 5. 

40. The compound of claim 10, wherein R is an electron pair, R' is H, W is O, and n is 6. 

41 . The compound of claim 1 , wherein the compound is selected from one of the following: 
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o- 




wherein, independently for each occurrence, 
R is H or CH2Y; 

R' is H, a sugar radical, or CH2Y; 
n is an integer from 1 to 100, inclusive; 
a is 1, 2, or 3; 
X is C(CH2Y)3; and 
Y is a protein binding group, 
wherein at least one Y is present in all compounds. 
42. The compound of claim 41 , wherein Y is a guanidinium ion. 
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43. A method of screening compounds for the property of inhibiting protein aggregation in 
solution, comprising: 

a) computing a set of parameters utilizing molecular modeling based on compounds 
known to have the property of inhibiting protein aggregation; 

b) applying those parameters to other compounds; and 

c) choosing the compounds that meet the criteria of those parameters. 

44. A method of preparing a compound having the property of protein aggregation inhibition 

in solution, comprising: 

a) computing a set of parameters utilizing molecular modeling based on compounds 
known to have the property of inhibiting protein aggregation; 

b) designing a compound having the property of protein aggregation inhibition in 
solution based on those parameters; and 

c) synthesizing the compound having the property of protein aggregation inhibition in 
solution. 

45. A method of classifying a compound as either inhibitory of protein aggregation in 
solution or not inhibitory of protein aggregation in solution, comprising: 

a) computing a set of parameters utilizing molecular modeling based on compounds 
known to have the property of inhibiting protein aggregation; 

b) applying those parameters to a compound; and 

c) classifying the compound that meet the criteria of those parameters as inhibitory of 
protein aggregation in solution. 

46. A method of determining the preferential binding coefficient, Fxp, of an additive in a 
protein solution, comprising: 

a) determining the phase space trajectories of the protein, solvent, and additive using 
molecular dynamics; 

b) calculating the distance, r, between the center of mass for both the solvent molecule 
and additive molecule to the protein's van der Waals surface; 
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c) determining the minimum distance, r*, at which no significant differences between the 
local (r = r*) and bulk density are observed; 

d) determining which molecules lie within the distance, r*, from the protein surface and 
classifying these molecules as the local domain; 

e) determining which molecules lie outside the distance, r*, from the protein surface and 

classifying these molecules as the bulk domain; 

f) determining the instantaneous preferential binding coefficient, FxpO), using the 
following formula: 



wherein: 

n"x = the number of additive molecules in the bulk domain; 
n'x = the number of additive molecules in the local domain; 
n"w = the number of solvent molecules in the bulk domain; and 
n'w = the number of solvent molecules in the local domain; and 

g) calculating the preferential binding coefficient, Fxp, as the time average of each of the 
values in step f) using the following formula: 



47. A method of suppressing or preventing aggregation of a protein in solution, comprising 
the step of combining in a solution a compound of the present invention and a protein. 

48. The method of claim 47, wherein the protein is a recombinant protein. 

49. The method of claim 47, wherein the protein is a recombinant antibody. 

50. The method of claim 47, wherein the protein is a recombinant human antibody. 

5 1 . The method of claim 47, wherein the protein is a recombinant human protein. 



rxp(t) = n"x - n'x (n"w / n'w) 




0 
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52. The method of claim 47, wherein the protein is recombinant human insuHn, 
recombinant human erythropoietin or a recombinant human interferon. 

53. The method of claim 47, wherein the solution is an aqueous solution. 

54. The method of claim 47, wherein the protein is a recombinant protein; and the solution 
is an aqueous solution. 

55. The method of claim 47, wherein the protein is a recombinant human protein; and the 
solution is an aqueous solution. 

56. A method of suppressing or preventing aggregation of a protein in solution, comprising 
the step of combining in a solution a compound of the present invention and a protein. 

57. The method of claim 56, wherein the protein is a recombinant protein. 

58. The method of claim 56, wherein the protein is a recombinant antibody. 

59. The method of claim 56, wherein the protein is a recombinant human antibody. 

60. The method of claim 56, wherein the protein is a recombinant mammalian protein. 

61 . The method of claim 56, wherein the protein is a recombinant human protein. 

62. The method of claim 56, wherein the protein is recombinant human insulin, recombinant 
human erythropoietin or a recombinant human interferon. 

63. The method of claim 56, wherein the solution is an aqueous solution. 

64. The method of claim 56, wherein the protein is a recombinant protein; and the solution is 
an aqueous solution. 

65. The method of claim 56, wherein the protein is a recombinant human antibody; and the 
solution is an aqueous solution. 

66. The method of claim 56, wherein the protein is a recombinant human protein; and the 
solution is an aqueous solution. 

67. A method of decreasing the toxicological risk associated with administering a protein to 
a mammal in need thereof, comprising the steps of adding to a first solution of a protein a 
compound of the present invention to give a second solution; and administering to a 
mammal in need thereof a therapeutic amount of said second solution. 
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68. The method of claim 67, wherein the protein is a recombinant protein. 

69. The method of claim 67, wherein the protein is a recombinant antibody. 

70. The method of claim 67, wherein the protein is a recombinant human antibody. In 
certain embodiments, the protein is a recombinant mammalian protein. 

7 1 . The method of claim 67, wherein the protein is a recombinant human protein. 

72. The method of claim 67, wherein the protein is recombinant human insulin, recombinant 
human erythropoietin or a recombinant himian interferon. 

73. The method of claim 67, wherein the first solution and the second solution are aqueous 
solutions. 

74. The method of claim 67, wherein the protein is a recombinant protein; and the first 
solution and the second solution are aqueous solutions. 

75. The method of claim 67, wherein the protein is a recombinant human antibody; and the 
first solution and the second solution are aqueous solutions. 

76. The method of claim 67, wherein the protein is a recombinant human protein; and the 
first solution and the second solution are aqueous solutions. 

77. A method of facilitating native folding of a recombinant protein in solution, comprising 
the step of combining in a solution a compound of the present invention and a 
recombinant protein. 

78. The method of claim 77, wherein the recombinant protein is a recombinant antibody. 

79. The method of claim 77, wherein the recombinant protein is a recombinant human 
antibody. 

80. The method of claim 77, wherein the recombinant protein is a recombinant mammalian 
protein. In certain embodiments, the recombinant protein is a recombinant human 

protein. 

81 . The method of claim 77, wherein the recombinant protein is recombinant human insulin, 
recombinant human erythropoietin or a recombinant human interferon. 

82. The method of claim 77, wherein the solution is an aqueous solution. 
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The method of claim 77, wherein the recombinant protein is a recombinant human 
antibody; and the solution is an aqueous solution. 

The method of claim 77, wherein the recombinant protein is a recombinant human 
protein; and the solution is an aqueous solution. 
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Figure 1 
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Figure 3 
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Figure 6 
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Figure 7 
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Figure 8 



0.1 
0.09 




1 

l\ 


Glycerol 

Urea 


0.08 








0.07 








0.06 




/I 




0.05 








0.04 




/ J \ 




0.03 




• •' 




0.02 


1 

/ 






0.01 


./ 






0 









-30 -20 -10 0 10 20 30 
Instantaneous Preferential Binding Coefficient 



-61 - 



Figure 9 
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Figure 10 
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